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Vortex- ant i vortex pairs in 2D easy-plane ferromagnets have characteristics of solitons in two 
dimensions. We investigate numerically and analytically the dynamics of such vortex pairs. In 
particular we simulate numerically the head-on collision of two pairs with different velocities for a 
wide range of the total linear momentum of the system. If the momentum difference of the two 
pairs is small, the vortices exchange partners, scatter at an angle depending on this difference, and 
form two new identical pairs. If it is large, the pairs pass through each other without losing their 
identity. We also study head-tail collisions. Two identical pairs moving in the same direction are 
bound into a moving quadrupole in which the two vortices as well as the two antivortices rotate 
around each other. We study the scattering processes also analytically in the frame of a collective 
variable theory, where the equations of motion for a system of four vortices constitute an integrable 
system. The features of the different collision scenarios are fully reproduced by the theory. We 
finally compare some aspects of the present soliton scattering with the corresponding situation in 
one dimension. 



I. INTRODUCTION 



The statics and dynamics of magnetic vortices is al- 
ready an old subject An increasing interest 
in the problem has arisen again §-§ which is con- 
nected with the synthesis and experimental study of 
new low-dimensional magnetic compounds such as two- 
dimensional magnetic lipid layers, organic intercalated 
quasi-2D layered magnets and HTSC-materials in the an- 
tiferromagnetic state. 

Magnetic vortices play an important role in easy-plane 
magnets. They are the main ingredients in the Kosterlitz- 
Thouless phase transition. At a finite temperature the 
density of vortices is large and they should give a consid- 
erable contribution to the dynamic correlations. 

We already have a detailed picture of the dynamics in a 
system with a small number of vortices. An isolated vor- 
tex in an infinite system can only move together with the 
background flux ||,^| . Two vortices interact and undergo 
Kelvin motion if they have opposite topological charges 
while they move one around the other when they have 
the same charge. We shall call the pair of two vortices 
with opposite topological charge a vortex-antivortex pair 
(V-A pair) . Such pairs can have the characteristics of a 
soliton in the sense that they move coherently with some 
constant velocity. Solitons of this kind have been nu- 
merically investigated in some magnetic systems PJlOf . 
Analogous V-A pairs have been studied in superfluids 
]TT| |l3| , nonlinear optics [Q and hydrodynamics E^j. 

In a system with a lot of vortices the picture becomes 
accordingly more complicated. If we suppose a dilute 
vortex gas, the average velocity is V ~ yfp ~ 1/L, where 
p is the density of vortices and L the average distance 



between them. The interaction energy is proportional 
to the logarithm of the average distance. This picture 
should be realistic above the Kosterlitz-Thouless temper- 
ature. At low temperatures vortex-antivortex pairs are 
expected to form. The energy of a pair is finite and it 
is proportional to the logarithm of its size. The interac- 
tion potential between them is inversely proportional to 
the second power of their size. One may be tempted to 
treat the vortex pairs as elementary weakly interacting 
particles. However, their dynamics is not Newtonian and 
most importantly they have an internal structure which 
may change during interaction. 

In a dense enough gas of vortices, interactions among 
traveling V-A pairs are unavoidable. In particular, any 
change in the number of vortices present in the system, 
or the number of vortices in equilibrium, should be a di- 
rect or indirect result of the scattering among V-A pairs. 
Our article is devoted to head-on and head-tail collisions 
between V-A pairs in easy-axis ferromagnets, i.e. to the 
case of zero total angular momentum of the system. Our 
study is both numerical and analytical. A collective co- 
ordinate theory is found to be particularly successful and 
provides the basis for a clear picture of the dynamics. 

Our results should also be relevant for a variety of other 
systems where V-A pairs have been found. Furthermore 
comparisons can be made to the well-studied soliton in- 
teractions in one space dimension. 

The outline of the rest of the paper is as follows. In 
Sec. II we give a short description of the system and an 
account of the dynamics of vortices and vortex-antivortex 
pairs. In Sec. Ill we present numerical simulations for 
collisions between V-A pairs. Sec. IV presents a theory 
which explains the features of the dynamical behavior of 
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vortex pairs. Our concluding remarks are contained in 
the last Sec. V. 



II. VORTICES AND VORTEX- ANTIVORTEX 
PAIRS 

We consider the classical two-dimensional Heisenberg 
ferromagnet with a uniaxial anisotropy of the easy-plane 
type. The corresponding Hamiltonian has the form 



H 



-J 



E 

(n,m) 



(i) 



where S„ denotes the spin variable at site n, S% is the 
third component of S n . The first summation runs over 
the nearest-neighbor pairs. The exchange constant J and 
the single-ion anisotropy constant fj are positive. We 
treat the spin S as a classical vector of constant length. 
Usually the magnetic anisotropy is small: /3/J^10~ 2 or 
1CP 3 . In this case we can use a continuum approxima- 
tion of the Hamiltonian. We define two fields to = S z 
and = axctan.{S v / S x ) , where S x , S y , S z are the carte- 
sian components of the spin. In terms of these variables 
the continuum version of the Hamiltonian reads 



n- fi 



dxdy 
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+ (1 - m 2 ) (V$) 2 + 



m 
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The coordinates x, y are measured in units of the "mag- 
netic length" Iq — J J I (3. to, $ are canonically conjugate 
fields and the equations of motion have the Hamiltonian 
form Q 



<9$ _ SH 

dt dm ' 
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Our numerical algorithm uses the formulation through 
the stereographic variable 



1 



1 + TO 

This satisfies the equation 
2H 
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where denotes the complex conjugate of f2. This vari- 
able was used in the solution of the Landau-Lifshitz equa- 
tion in one dimension since, in terms of it, the soliton 
solutions attain their simplest form |17|,p8|. 



In studying statics and dynamics for the above model, 
the topological density 
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dm d$ dm 9$ 
dx dy dy dx 



(7) 



is a most useful quantity. It has been called the "local 
vorticity" since it plays here a role analogous to the or- 
dinary vorticity in fluid dynamics |^ 19 1. The integrated 
topological density 



r 



7 dxdy 



(8) 



is an invariant and takes values which are integral mul- 
tiples of 2tt for the vortex solutions that we shall discuss 
here. 

The vorticity density enters the definitions of the linear 
and angular momentum of the theory which read Q 

P l i = ~£ t ii> x v jdxdy, fi,v = l,2, (9) 



{x 2 +y 2 ) 1 dxdy. (10) 



The role of the total vorticity in the dynamics can be 
also appreciated through the definition of the so-called 
" gyrocoupling vector" |20| 



G = -iY 



(11) 



where z is the unit vector in the third direction. The 
gyrocoupling vector enters the equations which describe 
the dynamics of vortices in a collective coordinate theory. 

We now turn our attention to the discussion of topo- 
logical excitations. We take as a boundary condition that 
the field <& is proportional to the polar angle <f> at spatial 
infinity. We then obtain vortex solutions which have the 
form @H1 



$ = n arctan 



y 



Y 



X ' 



= /(|r-R| 



(12) 



(13) 



where k — ±1,±2,... will be called the vortex number 
and (X, Y) is the position of the vortex center. In the 
following we call the vortices with k < antivortices. 
The magnetization field to for a vortex can be found nu- 
merically and has the following asymptotic behavior fl: 



to = A [1 — ar 2 ], r^O, 
to = A 6 exp(— r)/^/r, r ^> 1 



(14) 



where r is the distance from the vortex center, a, b are 
constants and A — m(r = 0) = ±1 we call the "polarity" 
of the vortex. The radius of the vortex is unity in our 
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units. We use in our numerical simulations only vortices 
and antivortices with k = ±1 and polarity A=l. 

The total vorticity of the vortices (|l2|,0|) is T = —2ttkX. 
The structure of the magnetic vortices ( ]12|jl3] ) is similar 
to that of the vortices in a non-ideal Bose gas jll| where 
the quantity (1 — m) is the density of the Bose particles. 
However, the ferromagnetic vortices differ in that they 
come with two possible values of the polarity A. 

The most impressive characteristic of the dynamics of 
an isolated vortex is that it is spontaneously pinned in 
an infinite medium. One can trace the reasons of this dy- 
namical behavior to their topological complexity which 
is reflected in the nonzero value of T ||. On the other 
hand, a vortex-antivortex pair undergoes Kelvin motion. 
This motion was studied in |2l],|2^] (see also j2j|) for a 
large vortex separation. In general, Kelvin motion sets 
in when the two vortices have opposite total vorticities: 
KlAl = -n 2 \ 2 . 

One can argue that the simplest topologically nontriv- 
ial objects which can be found in free translational mo- 
tion should have the form of a vortex-antivortex pair with 
r = 0. A conclusive numerical and analytical study in an 
easy-plane ferromagnet was given in JoJ where the profiles 
of coherently moving structures were numerically calcu- 
lated. There is a branch of solitons with velocities varying 
from zero to unity, which is the velocity of spin waves in 
the medium in our units. For small velocities the solitons 
have indeed the form of a V-A pair with a large separation 
L between the vortex and the antivortex. Their velocity 
is inversely proportional to the distance between them 
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(15) 



When the velocity approaches that of spin waves in the 
medium the vortex-antivortex character of the soliton is 
lost. The transition occurs at v ~ 0.78 in the sense that 
above this velocity the spin does not reach the north pole 
at any point. 

Similar ideas have a long history in hydrodynamics 
where V-A pairs have been studied theoretically and ex- 
perimentally in 2D flows. It is known that a V-A pair 
undergoes Kelvin motion in a direction perpendicular to 
the line connecting the centers of the vortex and the an- 
tivortex and the velocity is inversely proportional to the 
distance between them [15 lql . However, all the studies 



were made in the limit of point vortices, that is all mov- 
ing objects had a clear vortex-antivortex character. We 
shall go beyond this situation in the present paper. 



III. COLLISIONS OF VORTEX-ANTIVORTEX 
PAIRS: NUMERICAL SIMULATIONS 

The investigation of the interactions among the travel- 
ing V-A pairs comes as a natural next step after the well 
established theories of the previous section. In the follow- 
ing, we approach the subject through numerical simula- 
tions. We simulate collisions of pairs which are initially 



moving along the same line, say the horizontal x-axis. 
We are interested both in head-on collisions, where the 
pairs move initially in opposite directions and in head-tail 
collisions, that is when they move in the same direction. 

In the first set of our simulations, we confine ourselves 
to V-A pairs with small velocities so that the topological 
characteristics remain distinct and the vortices in the pair 
retain their identity. In particular, the distance between 
the vortex and antivortex centers is larger than the size 
of a single vortex (L > 1 in our conventions) . When the 
sizes of the simulated pairs are large, that is the distance 
between the vortex and the antivortex is large (L 1), 
we use as an initial condition the ansatz 



a = 



' i-fi 
i + fi 



(16) 



where and $i are the functions in Eqs. (|12j|13|). They 
represent vortices which are centered at Ri = (Xi,Yi) 
and have vortex numbers Ki. In particular we choose 
Ki = — k 2 , Xi=X 2 = Xi, Y 1 — —Y 2 = yi for the first pair 
which has a size L\ — 1y\ . We take accordingly ft 3 = — K4, 
and X3 — X4 = x 2l Y3 — — 14 = y 2 for the second pair 
which has a size L 2 = 2y 2 (we suppose 2/1,3/2 > 0). The 
above ansatz represents two V-A pairs which are moving 
on the x-axis and are set in a collision course, as long as 
the size of each pair is smaller than the distance between 
them: L%, £2 *C S = \xi — x-z\. As an alternative to the 
ansatz (|l6|) we also use the product ansatz of two V-A 
pair solitons of ref. M. This ansatz resembles in its gross 
features the ansatz (161). 

We perform the numerical simulations on a 500 x 500 
mesh with a uniform lattice spacing, typically h = 0.2. 
The time integration is performed by a fourth order 
Runge-Kutta routine. 

In PU one of the present authors has investigated the 
interaction process of two identical V-A pairs which col- 
lide head-on. It was found that during the collision the 
vortices exchange their partners and two new pairs are 
formed which are scattered at right angles. The trajecto- 
ries of the magnetic vortices are similar to those in colli- 
sions of V-A pairs in hydrodynamics (ljjl^l and they con- 
form to a good accuracy to the formula 1/Xf + \jYf — 
const, obtained in the 19th century p5| , pq |. 

Here, we consider collisions of two V-A pairs with dif- 
ferent velocities and consequently with different sizes. 
The results of our numerical simulations are summarized 
in Figs. [j],[||2] and ||. Fig. [l] presents head-on collisions and 
Fig. head-tail collisions through contour plots of the 
field m{x, y) at six snapshots during the collision process. 
Figs. |^ and ^ present the corresponding orbits of the in- 
dividual vortices during the interaction process. We have 
traced the center of every vortex which was considered to 
be at the point where the field m= 1. The result of each 
process depends essentially on the difference between the 
sizes of the two V-A pairs. For head-on collisions we use 
«i = K4, = 1, K2 — = — 1 and for head-tail collisions we 

USe Kl —K^ — 1, ft 2 = K4 = — 1. 
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(a) (b) (c) 

FIG. 1. Head-on collisions between vortex-antivortex pairs. We plot here contours of the third component of the spin using 
the levels 0.1,0.3,0.5,0.7,0.9. The numbers 1,4 denote the vortices and the 2,3 the antivortices. In (a) we have the case of a 
small difference between the linear momentum of the two pairs and the vortices exchange partners and scatter at an angle. 
In (b) the difference in momentum is larger. The pairs exchange partners, follow a looping orbit and finally rejoin the initial 
partners and travel along the initial direction of motion. In (c) the momentum difference is large. The two pairs pass through 
each other. Note that the boxes presented here have dimensions 40 x 40 while the simulations were done in a space 100 x 100. 
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Fig. 0a presents the head-on collision of two pairs with 
a small difference between their sizes. We use here the 
V-A pair solitons which have been numerically calculated 
in 1^ . The initial ansatz in our simulation is the product 
ansatz of two such pairs. We have taken the size of the 
left pair L\ ~4 which corresponds to a velocity v\ =0.27 
and the size of the pair on the right L2 — 6.3 which cor- 
responds to V2 = 0.15. The initial separation measured 
on the x-axis is (5 = 25. In the scattering process the vor- 
tices exchange their partners and two new identical pairs 
are formed which are scattered at an angle. Varying the 
sizes Li, L 2 we observe that the angle tends to 90° as the 
difference in the velocities (and momentums) of the two 
pairs is getting smaller. This proce ss g eneralizes the 90° 
scattering of two identical solitons |24| . 

According to (^), the larger soliton (pair on the right) 
has also a larger linear momentum. Conservation of the 
total momentum implies that each of the resulting iden- 
tical pairs has a non-zero cc-component of the momentum 
and a velocity to the left. Fig. 0a shows the trajectories 
followed by each vortex and antivortex. Open circles de- 
note the centers of the antivortices and filled circles those 
of the vortices. 

We defer for later the case of an intermediate differ- 
ence between the sizes of the two pairs and discuss first 
the case of a large difference (Fig. 0c and 0c). In the 
latter case one expects that the vortices which belong to 
different pairs will interact loosely with the vortices of 
the other pair. As a result, the two pairs are expected to 
travel almost undistracted. We use the ansatz (16) with 
the parameters L\ = 4, L2 = 1L\ = 28. The initial separa- 
tion of the pairs is 5 = 30. The result is close to expecta- 
tions, that is the small pair passes through the large one. 
The distortion in the trajectories should become smaller 
as the size of the large pair becomes larger. This case is 
thus analogous to soliton interaction in one-dimensional 
integrable systems, as has been noted by Aref PTj. 

In order to explore further this analogy we have plotted 
in Fig. 0a the x-coordinate of the two pairs as a function 
of time. (The data of Fig. 0c correspond to the data of 
Fig. [|a only until time=345). We observe that the fast 
pair experiences a delay during the interaction with the 
slow one, which results in a negative shift (in comparison 
to the free motion) in its ^-position after the collision. On 
the other hand, the slow pair is accelerated during the in- 
teraction and thus gains a positive shift in its x-position. 
The present observation should be contrasted to the sit- 
uation in one-dimensional soliton collisions where a pos- 
itive shift for both solitons is observed in all models. 

An intermediate situation between those in Figs. 0a 
and gc is presented in Fig. 0b. The parameters in the 
initial ansatz are L\ =4, L2 — 6L1 = 24, 6 — 34. The 
pairs initially exchange partners during the scattering 
process and form new V-A pairs just as in the simula- 
tion in Fig. 0a. However, after some excursion the new 
pairs approach each other again, exchange partners once 
more and the initial pairs re-emerge traveling along their 
initial direction of motion. The trajectories followed by 



the vortices are depicted in Fig. 0b. A similar scenario 
for point vortices in hydrodynamics has been discussed 
in (28|. 
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FIG. 2. The orbits of the vortices and antivortices of Fig. [I] 
during their head-on collision. The filled circles denote the 
position of vortices and the open circles the position of the 
antivortices at successive and equal times intervals. The num- 
bers 1,2,3,4 denote the initial position of the vortices and 
antivortices. 
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(a) (b) 

FIG. 3. Head-tail collisions between vortex-antivortex 
pairs. We plot here contours of the third component of the 
spin using the levels 0.1,0.3,0.5,0.7,0.9. The numbers 1,3 de- 
note the vortices and the 2, 4 the antivortices. In (a) the mo- 
mentum difference is large and the pairs pass through each 
other. In (b) we have a propagating quadrupole state. The 
boxes have dimensions 40 x 40, the simulations were done in 
a space 100 x 100. 
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FIG. 4. The orbits of the vortices and antivortices of Fig. [B| 
during their head-tail collision. The filled circles denote the 
position of vortices and the open circles the position of the 
antivortices at successive and equal times intervals. The num- 
bers 1,2,3,4 denote the initial position of the vortices and 
antivortices. 

In the next set of simulations we explore the situation 
of a head-tail collision. That is, both the slow and the fast 
pair move to the same direction (to the right in Figs. |^a, 
||a). The parameters here are L\ = 8, L 2 = 20, 5 = 25. 
The two pairs pass through each other, but some differ- 
ences to the case of Fig. |l|c should be pointed out. In 
Fig. ||b we give the ^-component of the trajectories of 
the pairs. We find that the fast pair is accelerated dur- 
ing the interaction (positive shift) while the slow pair is 
decelerated (negative shift). (See, however, the relevant 
remarks in the next section.) 

The last simulation, presented in Figs. ||b and |]b, in- 
cludes two identical V-A pairs traveling along the same 
direction. It can be considered as a limiting case to that 
of Fig. ||a when the sizes of the pairs are equal. The pa- 
rameter values are Li=8, £2 = 8, 5 = 6. The system can 
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be also viewed as a vortex- vortex pair and an antivortex- 
antivortex pair. Both pairs rotate while at the same time 
the magnetic quadrupole which is formed is propagating 
along the ir-axis. A similar "leap-frogging" motion was 



9| and has also been 
Since the magnetic 



studied in hydrodynamics by Love 
observed with two vortex rings [[50 
quadrupole is characterized by two parameters (the ve- 
locity and the internal frequency) it can be considered as 
an analog of a breather |27|^8| . We can make an estimate 
of the mean velocity of propagation of the quadrupole. 
Suppose that L is the mean distance of the pair of vor- 
tices from the pair of antivortices. If this is large com- 
pared to the distance between vortices of the same kind 
we may consider the quadrupole as two dipoles on top of 
each other. Then, a straightforward generalization of the 
results of M gives 



2 

I' 



(17) 



The rough estimate says that the quadrupole propa- 
gates with twice the velocity of a single V-A pair. In 
the present case, we have L — 8 which implies a ve- 
locity v ~ 0.25. Indeed, our simulations show that the 
quadrupole in Figs. ||b, ||b, propagates with a mean veloc- 
ity v q ~ 0.24 while a single V-A pair with a size L = Li=8 
has a mean velocity Vd — 0.125. 

The simulations of this section give an overview The 
richness of the above results shows that exploring numer- 
ically the different scenarios that occur during the inter- 
action of two V-A pairs is a rather cumbersome task. We 
note that in the processes that we have been studying all 
the vortices retain their identity during collision. This 
implies that an analytical calculation based on collective 
coordinates should be successful in reproducing the nu- 
merical results and should also provide an overview of 
the observed phenomena. This task will be taken up in 
the next section. 

In the remainder of the present section we continue 
with simulations of collisions of solitons which do not 
have a distinct vortex-antivortex pair character. We use 
here again the semitopological solitons found in ||. An 
example is given in Fig. pi In the initial state, in the 
first entry of the figure, the V-A pair on the right has a 
velocity v-i = 0.1 while the other soliton on the left has 
a large velocity, specifically here v\ = 0.9, and has no 
vortex-antivortex character. The clear numerical result, 
shown in the remaining five entries of the figure, is that 
the fast soliton is split at collision time in a vortex and 
an antivortex. The outcome is two identical V-A pairs 
which are scattered at an angle. 

We perform a series of simulations of collisions between 
solitons, where one of them has a definite velocity - we 
chose i>2 = 0.1 - while the velocity of the other varies in 
the range 0.1 < v\ < 1. The results are summarized in 
Fig. [7] which gives the cosine of the scattering angle 6 
as a function of the velocity v\ of the fast soliton (filled 
circles connected with a solid line in the figure). 
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FIG. 5. (a) shows the positions on the :r-axis of the two 
pairs (12) and (34) as a function of time, for the simulation 
of Figs. hlc. (b) shows the positions on the a>axis of the two 
pairs (12) and (34) as a function of time, for the simulation 
of Figs. Ha. See text for details. 

The scattering angle is measured (in Fig. ||) from the neg- 
ative horizontal axis. We have 6 = tt/2 when the velocities 
of the two solitons are equal, that is when v\ =V2 = 0.1. 
It is then monotonically decreasing as v\ increases until 
the value v\ ~0.9. 

When the velocity of the second soliton becomes 
greater than the value v\ ~0.91 we face a new scenario. 
The soliton is initially split into a V-A pair, then fol- 
lows a loop and finally it re-joins its initial partner. We 
eventually obtain the picture of two solitons which have 
passed through each other. The scenario resembles that 
of Fig. [ljb. The difference here is that, well after the 
collision time, the fast soliton is destabilized. When the 
velocity of the fast soliton is increased, the loop that 
the vortices follow during collision becomes smaller. For 

i>i ^0.97 we have no loop any more, and the fast soliton 
passes through the slow pair. 

The value «2 = 0.1 used in the simulations in Fig. |?j 
is not special. We can obtain results similar to those in 
Fig. using a different «2, e.g. i>2 = 0.2 or 0.3. 
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FIG. 7. The cosine of the scattering angle for a series of 
simulations where the slow soliton has a velocity V2 =0.1 
and the fast one 0.1 < vi < 1, is given by the dots which have 
been connected by a solid line. The dashed line results from 
conservation of energy and linear momentum. The dotted 
and the dot-dashed lines give the result of formula ( |35| ) for 
two different methods of calculating the soliton lengths (see 
text for details). 



IV. COLLISION OF VORTEX- ANTIVORTEX 
PAIRS: ANALYTICAL DESCRIPTION 



FIG. 6. Head-on collision between a slow V-A pair with 
velocity t>2 = 0.1 (right in the first entry) and a fast soliton 
with velocity v\ — 0.9 (left in the first entry). We present six 
snapshots at times t — 0, 22, 44, 66, 88, 132. The fast soliton 
is split in a vortex and an antivortex at the time of collision. 
The vortices exchange partners and two new pairs are formed 
which are scattered at an angle. 



The value of v± until which scattering at an angle occurs, 
seems to be close to the value «i =0.9 in these cases, too. 
For increasing vi the minimum scattering angle, that we 
can obtain by simulations, increases. 

The simulations of this section give an overview of 
the possible scattering scenarios. However, the picture 
will not be complete until we obtain an analytical under- 
standing. The theory which we present in the next sec- 
tion accounts for the simulations presented in Figs. [l],|3| 
and provides a reasonably satisfactory understanding. A 
step towards the understanding of the results of Figs, 
will also be taken. 



A full analytical investigation of interactions between 
vortex pairs appears to be quite complicated. Suffice it 
to say that no analytical formula is known for a single 
vortex pair soliton. However, one can employ a rigid 
shape approximation and suppose that each vortex is a 
coherently traveling structure and is also well separated 
from all the others. Then the dynamics of the system of 
vortices reduces to that of their centers R^. The latter 
obey the equations ]2C|fl 

^xG 1 = F,, (18) 

where Gi is the gyrocoupling vector (|ll]) for vortex i and 
F, = ~dE /9R,. The quantity E is the interaction en- 
ergy between the vortices, therefore we call Fj the force 
on vortex i exerted by the other vortices in the system. 

To make further progress we need the form of E, that 
is we need to know the field of the vortices. Under the 
assumption that the vortices are "quasi-static" the field 
of a single vortex is approximately given by Eq. ([l2],[l3]) . 
One should note that the energy is finite only when 

= 0, (19) 

i 

where Kj is the vortex number of a single vortex. Wc 
shall study only this situation here. Then the energy E 
of vortex interaction has the form pll] 
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E ~ -2tt kj ln(|Ri - Rj|), 



(20) 



where a constant has been omitted. 

We put R 4 = (Xi, Yi) in Eq. (|l|) and obtain 

Kj I 



"dT ~ ~ ^ Kj ' (R~ 



dYi _ \ (X^ — Xj) 
<ft~ ~ ^ Kj (R — R,) 2 • 



(21) 



These are the same as the equations of motion of point 
vortices in hydrodynamics |l5| , [l6|] when the hydrody- 
namic " vortex strengths" , which correspond to «j in the 
present system, are ±1. 
One can now see that 



i i 



as well as 



(22) 



(23) 



are conserved quantities. Eqs. ( |22| ) give the two compo- 
nents of the total momentum and Eq. (^3|) gives the total 
angular momentum of the system. They can be derived 
from the formulas fljlOl) when the rigid-shape approxi- 
mation is used. We further note that Eqs. ( |l8| ) can be 
derived as the Hamilton equations associated with the 
Hamiltonian ( po|) with the conjugate variables Xi and 

Pf = 2TTKiYi. 



In general, system (21) is integrable only for a system 
of three point vortices. However, for a system of two V-A 
pairs, for which condition (|l^) is satisfied, Eqs. ( |2l| ) can 
be integrated, when the pairs are initially moving along 
the same line [^2| . This is a fortunate situation since we 
shall be exclusively concerned with such systems in the 
present paper. 

The approximations employed so far seem to be quite 
crude. For instance we expect the rigid-shape approxi- 
mation to be valid when the distances between the vor- 
tices are much larger than the size of the out-of-plane 
structure of each vortex (unity in our units) . This could 
severely restrict the applicability of Eqs. (pi]). However, 
the numerical simulations, of the previous section, indi- 
cate that the results of the present approximate theory 
could be qualitative correct even for situations beyond 
the applicability limits of the theory. 

In the following we give the equations of motion of two 
V-A pairs. The vortex and antivortex of the first pair 
are placed at positions (xi,±yi) and have k = ±1 . We 
denote this pair schematically as (12). The second pair 
is at (x 2 , ±3/2)1 has K = ^fl, and is denoted (34). (cf. 



top entries of Fig. [j].) We choose the polarity A = 1 for 
all vortices. The system represents two V-A pairs in a 
head-on collision course. From Eqs. ( pl| ) we derive the 
four equations of motion 



dyi 
dt 



dx 2 
~~dT 



dm 

dt 



2/i - 2/2 



1 

2yx (xi - x 2 ) 2 + (2/1 - y 2 ) 2 

2/1 + 2/2 

(xi -x 2 ) 2 + (2/1 +2/2) 2 ' 

_ Xi — X2 

(xi - x 2 ) 2 + (2/1 - 2/ 2 ) 2 

xi - x 2 

(xi - x 2 ) 2 + (2/1 + 2/2) 2 ' 



1 



+ 



2/i ~ 2/2 



2y 2 (xi - x 2 ) 2 + (2/1 - 2/2) 2 

, 2/1+2/2 

(xi - x 2 ) 2 + (2/1 + 2/2) 2 

xi - x 2 

(xi - x 2 ) 2 + (2/1 - 2/2) 2 

xi - x 2 
(xi - x 2 ) 2 + (2/1 +2/2) 2 ' 



(24) 



(25) 



(26) 



(27) 



The system ( [24| - |27j) is completely integrable since there 
are four independent conserved quantities. We shall use 
the energy and the x-component of the linear momentum: 



2/12/2 



(xi - x 2 ) 2 + (2/1 - y 2 ) 2 _ _ (o) (o) 



(xi -x 2 ) 2 + (2/1 +2/2)" 



(0) (0) 

2/1-2/2 = 2/1-2/2 ■ 



2/i 2/2 



(28) 



(29) 



We suppose that the pairs are initially on the x-axis and 
at an infinite distance from each other. Then, yf*\ 
denote the 2/-coordinates of the vortices at time t = —00. 



We define L\ = 2y t f^ 1 and L 2 = 2y 2 v> as the sizes of the 
pairs. We take L 2 > L\ while L\, L 2 ^> 1. Therefore we 
have for the velocities: 



. = o..(°) 



1 



1 

V2 = T 2 



(30) 



that is, (12) is the "small" and fast pair and (34) is the 
"large" and slower pair. Eqs. (p4]-^7|) were studied within 
the framework of hydrodynamics in |33|,^8 29 1. It was 
shown that the behavior of the pairs during scattering 
depends on the ratio a — L 2 jL\ = v\/v 2 . There are the 
following three cases: 

(a) when 1 < a < a.\ = 3 + 2^2 — 5.83, the pairs change 
partners during the process and scatter at an angle 
(Fig. |l|a). 

(b) In the intermediate case, when ol\ < a < a 2 — 

(V2 + vV5- 1)/(V2- 1)~8.35, the vor- 

tices change partners during scattering but later 
they rejoin their initial partners and travel along 
the initial direction of motion (Fig. [jib). 



9 



(c) When a > ct2 the fast pair passes through the slow 
one (Fig. [l]c). 

In Fig. H we give a schematic representation of the differ- 
ent regions in the (vi , V2) plane where the three scattering 
processes occur. 

In case (a) the scattering process can be described by 
(12) + (34) -> (13) + (24) which is a schematic repre- 
sentation of the change of partners during collision. In 
the limit of two identical pairs one can use Love's |2^] 
equation (E8h to obtain the following well-known result 



vl 



(yi 0) ) 



(31) 



which is in good agreement with the data of the numer- 
ical simulation of the Landau-Lifshitz equation j24j] . In 
the general case of two different ingoing pairs, the outgo- 
ing pairs will actually be identical. Eqs. ( p8|]29| ) give for 
the size of the two identical outgoing pairs 



J out 



(32) 



This formula is known for point-like vortices in 2D hy- 
drodynamics of incompressible fluids pi . 

The angle of scattering can be found if we use the laws 
of conservation of energy and linear momentum: 



E = 



E, +Eo 



P cos 9 = 



Pa -Pi 



(33) 



where E is the energy and P the absolute value of the 
linear momentum on the a;- axis of each of the final pairs. 
The angle is measured from the direction of motion of 
the slow pair. We obtain 



cos ft = 



Pi 



2P 



(34) 



The calculation of P requires the use of the energy conser- 
vation law and the energy-momentum dispersion relation 
which is known numerically . 

Fig. shows the numerical results (filled circles con- 
nected by lines) for the cos 9 as a function of v\ , keeping 
a constant V2 = 0.1. We have v\ > v%. The dashed line 
results from Eq. (|34|) and it is a very good approximation 
of the simulation results until V2 ~ 0.78. At this value 
the dashed line has a maximum. Until this critical value 
of i>2 the collision is elastic in the sense that almost no 
energy is dissipated in radiation. We also have to men- 
tion that the approximation shown by the dashed line is 
getting worse as the velocity v\ is increasing. 




FIG. 8. A schematic overview of the different possibilities 
for head-on and head-tail collision between two soliton pairs. 
vi and V2 are the velocities of the two pairs. The regions a,b 
and c correspond to to the cases a,b and c of Fig. |l| The 
region d corresponds to the case of Fig. j^a and the scenario 
of Fig. ^b occurs for velocities on the line OE. 



and p^) we have 



In the approximation of Eqs. (|2 
for the scattering angle: 



1 



cos ft = 



(35) 



The limiting cases are 8 = ir/2 for a= 1 and 9 = tt when 
a = a.\. Some care is needed in interpreting the angle 
ft = 7r in the last limiting case. In this case the new pairs 
which emerge follow the parabolic orbit 



Xl 



2y? 



vl 



2(l + x/2)yi' 



(0)' 



(36) 



The above results, should be a good approximation in 
the limit of large Pi, L>2- 



There is a simple way to apply Eq.( 35) when the two 
colliding V-A pairs are slow and the vortex and antivor- 
tex in a pair are well separated. The size of each pair can 
be taken to be the distance between the two points where 
the spin variable reaches the north pole and its velocity 
is the inverse of it. This method gives fairly good results 
for simulations with pairs of the size used in Fig. [j]. 

A comparison of Eq. ( |35|) with numerical results is 
given in Fig. ^. The cosine of the angle 6 for 1 < a < ati, 
is plotted by a dotted line. Fig. fj] shows cos 9 as a func- 
tion of vi, for V2 = 0.1. (not as a function of a). We 
have vi > v 2 and we consider a = Vi/v-2- The dotted 
line is a poor approximation to the simulation results for 
a > 4 (v2 > 0.4). The deviation from the numerical 

points is even qualitatively wrong already for V2 ~ 0.583. 

We recall at this point that in our simulations of the 
previous section we obtained results which could be un- 
derstood as interaction processes between two dipoles 
even in the cases where the V-A pairs have no apparent 
dipole character. Exploiting this remark we assume that 
the V-A pair solitons (even those with large velocities 
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and no apparent vortex-antivortex character) are dipoles 
with a length L and a charge q. To be sure, in the limit 
that the velocity goes to zero, L should go to the simple 
definition of length described in the previous paragraphs 
and q should go to the vortex number k = ±1. In gen- 
eral, however, we use a generalization of Eqs. (|2^) and 
( 15|) and write for the momentum and velocity of such a 
dipole: 



P = 2-KqL, v=y. (37) 
The length and charge of the dipole are then given by 



L = 



P 

2nv ' 



q = 



(38) 



where the values for P, v can now be taken from the Table 
of ref . . In the limit v — > we use the relation Pv = 2n 
to find that L = 1 jv — ► oo and q = ± 1 , which is the ex- 
pected limit. In the opposite limit v— *1, an asymptotic 
analysis || gives that P~ 1/ Vl — v 2 — > oo, that is an in- 
finite length of the dipole and also q ~ l/Vl — v 2 — ► oo. 
For intermediate velocities the length L is finite and it 
reaches a minimum for v ~ 0.87. The values for the 
charge q are relatively close to unity for v < 0.5 and they 
increase rapidly for v > 0.9. 

The interesting result is presented in Fig. 0by the dot- 
dashed line. It is obtained by applying Eq. (|35|), where 
for a we substitute the ratio of lengths of the two V-A 
pairs given from Eq. (|3^). The curve is compared quite 
good with the results of the simulations even for large 
velocities. However, after v 2 = 0.9 the present approach 
fails completely. Indeed, we have no reason to believe 
that an approximation based on dipoles would be cor- 
rect in the limit v — > 1. 

We now turn to case (c) where the initial pairs survive 
throughout the process. For a large enough an expansion 
gives, with an error 0(1/ a 4 ), the trajectories 



X\ = Vi t ■ 



arctan (2v\V 2 t) 



x 2 



-v 2 t 



1 2 {2viV 2 t) 
a 2 v 2 1 + (2viv 2 t) 2 
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— arctan(2t>iV2i) 
or v 2 



yi = yi 



(0) 



V2 = y 2 



(0) 



1 



a 1 + (2v!V2ty 



1 



1 + 



a 2 1 + (2viv 2 ty 



(39) 



where v\ , v 2 are given in Eq. ( p0[ ) and the dots in the 
second equation stand for a lengthy term of order \jo? 
which we shall not need in our analysis. It follows that, 
in this approximation, the fast pair (12) deviates strongly 
from the rectilinear motion. The shifts in the positions 



00 to 00, are 
3 for 



of the vortex pairs as time varies from 
Axi = — An/vi for the fast pair and A^2 = 47r/w2a 
the slow one. (The terms in the dots do not give any 
shift in the pair position when times goes from —00 to 
00.) The shifts are measured relative to the direction of 
motion of each pair. The results of the simulations for 
the vortex system show that the fast pair is decelerated 
during the collision as if it is repelled by the slow one, in 
agreement with Eq. (|39|). The simulations further show 
that the slow pair is accelerated during the collision as if 
it is attracted by the fast one. The final result is a pos- 
itive and negative shift in the positions of the fast and 
slow pair, respectively, a phenomenon which has not been 
observed in head-on soliton collisions in one dimension. 
The obtained trajectories are similar to those for the 2D 
Euler equation B4|. 

In case (b), a± < a < ct2, the scattering process is rep- 
resented by the scheme (12) + (34) -> (13) + (24) -> 
(12) + (34). This means that the vortices exchange part- 
ners at a first stage of the scattering, then the new pairs 
follow a looping orbit and at the final stage the initial 
partners rejoin and travel along the initial direction of 
motion (cf. Figs. |l|c, ||c). The loop becomes the parabola 
(|36| ) in the limit a. = ct\ and it is a cusp when a^a 2 - 

We now turn to the head-tail collision. It corresponds 
to the area denoted (d) in Fig. || and the relevant numer- 
ical simulations have been given in Figs. |^ and |]. The 
case has been studied within a hydrodynamical context 
As Love showed, only a slip-through mo 
) and a leap-frogging motion (Fig. ||b) can 



in |2^_ 
tion (Fig. 
occur. 

Suppose that vortices 1 and 3 have k — 1 while 2 and 4 
are antivortices and have k = — 1 (cf. Fig. ||a). The equa- 
tions of motion arc modified as follows: the left hand side 
of (Si], ^5]) and the first term on the right hand side o f 
( p6| ) change their signs. The conserved quantities 
now read 




(xi - x 2 ) 2 + (yi +y 2 ) 2 
(xi - x 2 ) 2 + (yi - 2/2) 5 



(yi 



t (0) 
(0) (0) (x 1 

Vi V* TTof 



(0) 



4 0) ) 2 



(0) 



(°) 

2/1 + 2/2 = y\ ^ 

M 



(2/i 



(0) 

•2/2 



2/f) 



±y<°>) and (xf\, 



(40) 



(41) 
±2/f). 



The pairs are initially at (x\ 
Using the above one can find that two pairs which start 
infinitely far apart will pass through each other for any 

value of yf^ < ■ 

For a large difference in the size of the pairs (a ^> 1) 
we find, by an expansion, the solution of the equations 
of motion (with an error O (-tt)): 



x\ = v\ t H arctan (2viV2 1) 

Vl 



1 2 (2viv 2 t) 

x 2 =v 2 t — 

a 2 v 2 1 + (2v 1 v 2 ty 
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1 4 

a 3 v 2 



arctan(2vii>2t) 



Vx 



1J2 



(0) 



(0) 
2/2 



1 - 



a 1 + {2viv 2 1) 2 



(42) 



The dots stand for terms of order O (^) which do not 
contribute to the shift in the position of the pair when 
times varies from — oo to oo. In the present case, the 
shift of the fast pair is to the direction of its motion 
Axi ~ Att/vi. The shift of the slow pair is also positive 
but small Ax 2 — 4:Tt/(v 2 cx 3 ). However, the numerical 
simulations of the previous section (Fig. ||b) have given 
a small shift for the slow pair opposite to its direction 
of motion. We believe that this should be a consequence 
of the second term in the second of Eqs. (^) which is 
dominant because of our small space and short integra- 
tion time. We note here again a difference of the present 
system with the situation in ID. For instance, in KdV 
the fast soliton acquires a positive shift and the slow one 
a negative shift due to a head-tail collision. 

A difference between the head-on and head-tail colli- 
sion for the case of large a, is that during the collision 
the size of the small pair increases in the head-on colli- 
sion case (Eq. (|3|) and Fig. Jl]c) but it decreases in the 
head-tail collision case (Eq. ( |42"| ) and Fig. 

In the limit of a small difference of the sizes of the pairs 
(a~ 1) the relative shift of the pairs grows and tends to 
infinity for identical pairs like Ax± — Ax 2 ~ 4/(ui— v 2 ). In 
the limit of two identical pairs the asymptotes of the solu- 
tion after the scattering have the form x 2 — Vt— \fi , x% = 
Vt + \ft and, at large distances, the distance between 
pairs goes as X\ — x 2 = 2\/t. It is interesting to compare 
this shift with that for solitons in ID systems where the 
shift is proportional to the logarithm of the difference of 
the velocities of the two solitons (Ax% ~rn(ui — U2)) and 
the distance between two identical solitons after collision 
goes like x\ — x 2 ^\nt. 

Two V-A pairs which have a finite distance between 
them, will pass through each other when the following 
relation holds 



(V. 



(o) 



(0)\2 

Vi ) 



(OK 2 



2(<4 0) -2/i 0) ) 



(oh 
Vi ) 



> (4 0) 



(43) 



In the opposite case, the pairs form a translating bound 
quadrupole state. The translation is accompanied by 
a rotation of the two vortices and the two antivortices 
around each other. This leap-frogging motion was first 
analyzed in ^9|. An example is given in Figs. ||b and [|b. 
The quadrupole state has lately acquired special interest 
due to its relation to breather modes pS} . 

We shall analyze some characteristics of the leap- 
frogging motion in the case that the two rotating vortices 



are well separated from the two rotating antivortices. We 
call L the distance between the two pairs. The two vor- 
tices rotate clock-wise and the two antivortices counter 
clock-wise. The translational motion of the quadrupole 
has the velocity given in Eq. (|T7|). With the further as- 
sumption that the distance 6 between the vortices of the 
same pair is large compared to the size of a single vortex 
(but still small compared to the size of the quadrupole: 
L,8 3> 1,L ^ (5), we find from Eqs. ( |39| ) the frequency 
of rotation 



2 



2 

l2' 



(44) 



V. CONCLUSIONS 

We have presented a detailed study for the head-on and 
head-tail collisions between vortex-antivortex pairs in 2D 
easy-plane ferromagnets. The V-A pairs are the simplest 
units which can be found in free translational motion 
in our model. However, they have a nontrivial internal 
structure which is responsible for their unusual behavior 
during interaction. Our study combines numerical sim- 
ulations, which yield a variety of interesting scattering 
scenarios, with a collective variable theory, which leads 
to an understanding of the main features of the scatter- 
ing scenarios. The change of partners between V-A pairs 
during interaction is the most remarkable effect in our 
case. It is due to the internal structure of the V-A pairs 
which has been fully taken into account. 

In a study of a system of a lot of vortices the mecha- 
nisms which have been described here should play a dom- 
inant role in the way to the final steady state. The re- 
laxation process should include a multitude of collisions 
of vortices and the formation of steady structures. 

From the point of view of soliton theory in two dimen- 
sions, our results are novel and can be compared with a 
variety of studies in the field. Some comparison is also 
done with soliton theory in one space dimension. 

The scattering of coherently traveling objects in two 
dimensions appears to be of interest in a variety of phys- 
ical systems. Objects similar to the V-A pairs studied 
here exist in systems in different fields of physics [[l2]-[l5) . 
Even in nonequilibrium systems ]35|,|36| studies similar to 
the present one have been performed and some similar 
results have been derived. 
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